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COMMENTS ON "CONTACT DIFFUSION INTERACTION 


OF MATERIALS WITH CLADDING" 
by James F. Morris 
Lewis Research Center 

i 

SUMMARY 

A Russian paper by A. A. Babad-Zakhryapina contributes much to the understanding 
of fuel, clad interactions and thus to nuclear thermionic technology. In that publication 
the basic diffusion expression is a simple one. A more general but complicated equa- 
tion for this mass transport results from the present work. With appropriate assump- 
tions, however, the new relation reduces to Babad-Zakhryapina’s version. 


A REACTOR MASS -TRANS PORT PROBLEM 

A major complication in nuclear power generators is the intermingling and reaction 
of the fuel with its immediate container. When the restraining wall is an essential part 
of the energy converter itself, this physicochemical interchange affects the output and life 
of the system directly. Then it becomes a crucial problem - as it is for in-core ther- 
mionics. 

Reference 1 focuses analytically and experimentally on part of this technical area, 
the "contact diffusion interaction of materials with cladding. " Through close theoretic 
and empiric coordination of its results, that publication adds substantially to the knowl- 
edge of mass transport in the reactor. In the first two sections Babad-Zakhryapina 
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develops the model and a simple diffusion equation to serve as the basis for further 

analyses. That initial presentation appears in this report in the appendix to assist the 

* 

reader and to obviate repeating the detailed introduction. 

The present derivation begins with the assumptions of reference 1 and produces a 
less approximate expression for the migration of fuel components through clads. Then 
a comparison of the new and the old descriptions indicates the applicability of Babad- 
Zakhryapina's version. 


A MORE GENERAL DIFFUSION EQUATION 


The basic fuel, clad diffusion problem is that outlined by Babad- Zakhryapina and 
quoted in the appendix. The solution sought is one for "low values of k and t, " as 
well as for very small x < l = kt. 

But before the new derivation a comment on the referenced one seems in order: 
The integral version of equation (A3) applies only for the location x = l, where 
dx/dt = dZ/dt = k and 3c/3x = 3c/aZ. Here c(x,t) = c^ = c(Z,t) = c(kt,t) = c(t) x _^. And 
the fully developed differential expression is 


V 


D dc ^x=Z 
k dt 


D 

k 



= c A k 


( 1 ) 


So equation (A 4) really becomes 


C A 




( 2 ) 


It describes the composition at the surface of the growing layer only - not the space, 
time variations of concentrations throughout the lamina. These effects are the objects 
of the following development. 

First from 


— = D V 2 c (3) 

at 


by Laplace transformation comes 


sc(x,s) = D^- c - (x ’ s) 
2 


3x 


(4) 


Then, 


c(x, s) = c(o, s)exp 



(5) 


Obviously c(o,t) is c(o,o)f(t), but f(t) hides effectively behind indefinite boundary con- 
ditions. 
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However, because k and t are slight, yet large compared with x, a "small 
system" approximation is possible (ref. 2). There, near homogeneity holds allowing 
separation of the variables: 


c(c,t) « X(x)T(t) 


( 6 ) 


Then Fick's second law becomes 


1 dT 1 


^d 2 x\ 


DT dt X ldx 2 ‘ 


(7) 


But this can hold only if each side equals a constant, -A . Thus, 


T = T q exp(-A^ Dt) 


( 8 ) 


results with 


X = A' sin Xx + B' cos Xx 


( 9 ) 


So the particular solution is 


c(x,t) « (A sin Xx + B cos Ax)exp(-A^ Dt) 


( 10 ) 


and the general expression is 


c(x,t) « 



01 ) 


For an approximate thin-film approach, though, the simplified solution suffices. 

As k - 0, t - 0, and x ^ kt prevail, sin Xx - 0, cos Xx - 1, and exp(-A Dt) - 1 follow. 
Therefore, the particular solution (10) reduces to 

O 

c(x,t) w c(o, o)cos(Xx)exp(-X Dt) (12) 


This shows that the corrected form (2) of Babad-Zakhryapina's relation implies con- 
centrations in the lamina because cos(Xx) — 1 typifies the small system version. Then 
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2 2 2 

a comparison of exponentials reveals the constant A must be k /D . And finally the 
particular solution becomes 


c(x,t) 


W 

c(o,o)cos — exp 
D 



(13) 


The variation of internal composition with distance, which the approximation of 
reference 1 lacks, appears as a cosine factor. Furthermore, the proper version of 
Babad-Zakhryapina’s expression (2) describes the changes in concentration with time at 
the interface (x = 0) in the small system answer (13): 


I 


c(o,t) 


c(o, o)exp 



(14) 


In turn, that provides a good estimate of the missing function in the more general 
Laplace development (5). With the transform of (14), 


c(o,s) 


c(o,o) 
(s + k 2 D) 


05) 


equation (5) yields 


c(x,s) » c(o,o)expi^) 

(s + k 2 /D) 


(16) 


Reference 3 indicates the inversion as 


c(x,t)~ c(Q ’ o) e- (k /D)t j 


e -i(k/ D )x erfc / x . lk r 

tNeUVDjx / x +ik fi\ 

^ c(o, 0) p -(k 2 /D)t 1 

y2v/bt y : 

d J V D / 

2 \ 


, k 
cos — x 


■- i sin k xV 1 - -1 J(-X- - ik t /l) - fcZgl^.. ik V^Pl 3 + Q ^VPt) - ikt/VD ? . j\ ( 17 ) 

D /l |Al\2|/bt f D J 1!3 2! 5 IJ 


+ |cos — 


k A/, 2 1/ X 


x + i sin — x)/l - 


VM\2 Dt 
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The real terms of equation (17) approximate Babad-Zakhryapina's diffusion model. 
Further reductions of (17) result from the restrictions of the lead paragraph of this 
section: "For Tow values of k and t, ’ as well as for very small x < 7. = kt, " 
sin kx/D and all terms except the first in the erfc expansions approach zero. Thus, the 
simplified form of the diffusion equation (17) emerges as 


c(x,t) 



(18) 


which verifies the small system result (13) and reinforces the corrected expression for 
Babad-Zakhryapina's approximation (2). > 

In fact, for vanishingly thin films x — kt — 0, cos(kx/D) — cos(k t/D) — 1, and 

O 

exp(-k t/D) — 1 all hold. The strengths of the exponential and the cosine effects, 
however, differ widely: 


d exp(-k 2 t/P) = _ k 2 /-k 2 t\ _ k 2 / k 2 t | k 4 t 2 \ _ k 2 

dt D \ D / D \ D d 2 ' ' / D 

d cos(k 2 t/P) _ k 2 G . n /k 2 t\ k 2 / k 2 t [ k 6 t 3 * * k 10 t 5 ( \ _ lA 

dt " D \ d )~ D\D 3!D 3 5 , d 5 ' ' d 2 


(19) 


( 20 ) 


2 2 

Thus, as kt - 0, exp(-kt /D) changes at least D/(k t) times more rapidly than 

cos(kx/D). This means, of course, that the exponential factor dominates the basic dif- 
fusion equation at "low values of k and t. " However, Babab-Zakhryapina's presenta- 

tion failed to show and explain this fact - or even to obtain the proper exponential de- 

pendency. 

One final point deserves mentioning: Although (13), (17), and (18) are adequate 

approximations for k — 0 and t - 0, substituting [l + cos(kx/D) - cos(k 2 t/D)] for 

cos(kx/D) in the equations is more satisfying. And, of course, the added quantity 
[l - cos(k 2 t/D)] — 0 is a very small time effect compared with exp(-k 2 t/D) as (19) and 

(20) testify. Its addition or deletion is well within the allowances of the present estima- 
tion. But this change makes (13) and (18) 


c(x,t) 


c(o,o) 1 


+ cos 




( 21 ) 
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Thus, at the interface c(o,t) = c(o,o)[2 - cos(k^t/D)jexp(-k^t/D), which is smaller than 

c(o,o) because the exponential decrease overwhelms unity minus the cosine. On the 

o 

outer surface (x = l = kt), however, c(?, t) = c(o, o)exp(-kZ/D) = c(o,o)exp(-k t/D). This 
is the correct result indicated by (1) and (2). So, even though (18) and (21) are both 
reasonable derivations from the model, the latter is better. 


IMPLICATIONS 

Thd" mass-transport equations developed here reveal space, time trends for fuel- 

atom concentrations in clads following the model of reference 1. These variations 

( 

given by 'a general but complicated relation (eq. (17)) correspond to those obtained with 
a small-system approximation. Furthermore, a final trival refinement yields an esti- 
mate (eq. (21)) that satisfies all boundary conditions - including the absolute solution for 
the outer-surface composition. 

Compared with the present findings, Babad-Zakhryapina's fuel- component distribu- 
tion predicts high local concentrations within the clad. Since his equation actually de- 
scribes only the outside condition (eqs.(A4), (1), and (2)), though, this deviation is in- 
tuitively obvious. But when the model operates at the extremes of its assumptions, the 
simple mass- transport expression of reference 1 is an admissible approximation. This 
is gratifying because Babad-Zakhryapina’s presentation represents a valuable contri- 
bution by eight Russian scientists on the "contact diffusion interaction of materials with 
cladding. ’’ 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 16, 1971, 

112-27. 
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APPENDIX - EXCERPT FROM REFERENCE 1, "CONTACT DIFFUSION 
INTERACTION OF MATERIALS WITH CLADDING" 
by A. A. Babad-Zakhryapina 

This published work discusses problems arising from diffusion interaction between materials 
and their claddings and the film of condensate formed on their surfaces. Since condensates 
practically always form on heated materials, the commencement of interaction is observed during 
the process of formation of the film. In this connection, the investigation of the characteristics 
of the diffusion interaction of the film with the materials should be divided into two stages: inter- 
action during the process of the formation of the film and the interaction during process of ex- 
tended isothermal heating. Detailed results of two stages are presented in this work. 

A characteristic of the first stage is the simultaneousness of the course of the two processes: 
condensation and diffusional transfer of the cladding material into a growing layer. The continuity 
of the condensate growth, as well as the law governing the growth of the layer itself, gives rise 
to the appearance of certain specific characteristics in the character of diffusion interaction 
(Section 1). 

The interaction of the film with the material during isothermal heating may occur at the ex- 
pense of a decrease of the content of one of the components of the cladding material and not as a 
result of the interchange reaction. It is borne out that geometric correlations have a substantial 
influence on the character of diffusion process in such a system. This is shown by the example of 
contact diffusion of monocarbide of uranium with molybdenum and tungsten (Section II). 

Besides these general questions concerning characteristics of diffusion interaction of deposited 
materials, investigation is also made of the correlation of the diffusional mobility of atoms with 
thermodynamic characteristics of materials. 

In conclusion, a study is made of data concerning the diffusion of carbon in zirconium carbide 
(Section III). 

1. CHARACTERISTICS OF DIFFUSION PROCESSES IN CONDENSATES, 

FORMED ON A HOT CLADDING 

by A. A. Babad-Zakhryapina and L. Gert 


1 • Introduction 

To create a well bonded layer of condensate with a cladding surface, it is common practice to 
utilize heating of the cladding during the process of condensation. During the course of the whole 
process, there occurs an inevitable diffusion interaction between the cladding material and the' 
condensate material, a characteristic of which is a diffusional penetration of the cladding material 
into the layer of condensate. On the other hand, the penetration of the condensate material into the 
cladding also, to a certain degree, has to be determined by the conditions of condensation. 

The diffusion interaction between the cladding and the condensate may lead to a change in the 
conditions of condensation, i.e. , a change in the character of dependence between the thickness 
of the layer and the time it takes to form, as well as have a substantial influence on the phase 


7 



formation in the condensate layer. Regardless of the fact whether these phenomena are desirable 
or undersirable in the process of formation of the condensate layer, they have to be taken into 
account and a method, based on prior knowledge, has to be developed by which the phenomena can 
be evaluated quantitatively. 

2. Unilateral Interaction of the Cladding Material with the Condensate 
Layer During the Process of the Layer Formation 


For the solution of this problem, it is necessary to know the growth law of the condensate 
layer. Subsequently, two laws will be examined: linear, i.e. , the growth of the layer with a con- 
stant speed and parabolic - with a variable velocity. 


The simplest case, which is encountered in the investigation of the diffusion interaction of the 
cladding with the condensate, is the unilaterial diffusion in the condensate layer. This case may 
be formulated in the following manner: on a cladding of material A at a temperature T, sub- 
stance B is placed according to the law x = f(t) (t - time); it is necessary to determine the con- 
centration of substance A (i.e. , c A ) on the surface of the layer of condensate. 


Here it will be assumed that the coefficient of diffusion A in B (D AB ) is considerably greater 
than the coefficient of diffusion of B in A (Dg^). For simplicity and ease of notation, let us con- 
sider that D t 


BA 


= 0 and 


°AB = D ‘ 


Boundary conditions for this problem have the form (fig. 1) on the interface cladding conden- 
sate. 



cM x=0 = c o 


on the surface of the layer 


c(x, t) x=£ = c A , where £ = f(t) 


From the general solution of the Fick equation for a unilateral source and above boundary condi- 
tions for a case when £ = 2b^Dt, we obtain: 


c = 


1 + F(b) 

where F(b) = b exp b 1 2 ^(b), 'F(b) - function of Gauss errors. 
Under linear law of growth of layer, when £ = kt 


(Al) 


c. = 


A 1 + F(k, t) 


(A2) 
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where 


The expression, determining under linear law of layer growth, is an approximation and 
may be utilized for calculations at low values of k and t. An accurate solution of this problem 
may be obtained by the solution of the equation of material balance on the boundary of cladding 
condensate, where 


-D — 
0x 


x=£ 


= c A k 


(A3) 


bearing in mind that d£/dt = k. 

Let us perform a substitution of variable k = £ . Such a substitution is possible provided that 


0c| _ dc A 
0t 1^ dt 


We will prove that the latter equivalence is valid only in the case if £ = kt. Differentiating equa- 
tion (A3) according to £ and comparing the obtained result with the Fick equations for point £ , 
we obtain 



As c A ^ 0 according to the problem, then (d/d£) (d£/dt) = 0 which gives £ = kt. 

Bearing in mind that c(£,t) = c A and -D(8c A /a£) = c A k, after transfer to variable "k", we 
obtain 


c(x,t) = c 0 



(A4) 


Let us prove that expression (A4) satisfies Fick's equation and boundary conditions. For this pur- 
pose, it is expedient to go to another system of coordinates X = £ - x (fig. 2). In the new system 
of coordinates, expression (A4) will be written as 


c(X,t) = c Q exp 



(A 5) 


By direct substitution, it is easy to be convinced that expression (A5) satisfies both the Fick 
equation and boundary conditions. 


From equation (A4), it follows that the distribution of the concentration in the layer of the con- 
densate, growing according to linear law, is pseudo -stationary. The form of the concentration 
curve in this case is determined only by the relationship k/D. 


q 



With the increase of the growth rate, the concentration of layers at any one given point, and 
consequently on the surface, decreases. The same conclusion is arrived at by analysis of approxi- 
mated correlation (A2). 

In any case of parabolic law of growth, c^ remains a constant independent of time while its 
value is determined only by parameter "b. " The above obtained correlations (A3) and (A4) may be 
utilized for determination of coefficients of diffusion. The most suitable for this purpose is cor- 
relation (A4). Conducting a linear interpolation of the concentration curve, which is permissible 
for lower values of kx/D, we obtain c(x, t) = c Q D/(D + kx). In this relationship, magnitudes of 
"k" and "x" are easily determined experimentally. Thus, for computation of D, it is only 
necessary to know the concentration at any one given point along the thickness of the layer. In the 
case where, as a result of diffusion, solid solutions are formed and if the concentration dependence 
of the lattice parameter is known, then the concentration on the surface of the layer can be deter- 
mined by X-rays. 
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Figure 1. - Diagram of the distribution of concentration in a layer of conden- 
sate in a case of unilateral single phase diffusion. (From ref. 1.) 


c 



Figure 2. - The connection between systems of coordinates showing the 
layer on the growth rate under the linear law for the growth of the 
condensate. (From ref. 1.) 
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